Collective polarization exchanges in collisions of photon clouds. 
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The one-loop "vacuum" Heisenberg-Euler coupling of four electromagnetic fields can lead to in- 
teresting collective effects in the collision of two photon clouds, on a time scale orders of magnitude 
faster than one estimates from the cross-section and density. We estimate the characteristic time 
for macroscopic transformation of positive to negative helicity in clouds that are initially totally 
polarized and for depolarization of a polarized beam traversing an unpolarized cloud. 
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Some non-linear aspects of vacuum electrodynamics 
have been tested in experiments on Delbruck scattering 
PJ , i. e. the scattering of a photon off of the Coulomb 
field of a nucleus, and in photon splitting 2], also in the 
nuclear Coulomb field. 

Essentially, these effects hinge on the one-loop effective 
Lagrangian density for processes in which four or more 
electromagnetic fields, of long wavelength compared to 
the electron Compton wavelength, come together,as de- 
scribed by the Heisenberg-Euler interaction Q 0], the 
fourth order term of which is, 



Lj= H 3 *-^[(E 2 -B 2 ) 2 + 7(E.B) 2 ], (1) 
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where a is the fine structure constant and m is the mass 
of the electron. 1 

The validity of the effective interaction term JJ), for 
long wave-length fields, can hardly be doubted. Nonethe- 
less, its confirmation in an actual photon-photon scatter- 
ing experiment would be a milestone of a kind. If one 
puts in the numbers for photon-photon scattering itself, 
the cross-section is far too small to be measured with 
current technology. Indeed the "light by light" scattering 
discussed in the very interesting experiment reported in 
ref. 5] was the reaction 7+7 — > e + + e~, and does not test 
vacuum QED at the one loop level. However, Kotkin and 
Serbo 6] have pointed out that a photon of one plane po- 
larization passing through a cloud of photons that are po- 
larized in a different direction in a frame in which the col- 
lisions arc head-on, will experience a polarization preces- 
sion with an angular frequency, T p = 4cv 2 n 7 o;aj c /(I5TO 4 ), 
where oj,uj c are the respective frequencies of the imping- 
ing photon and the cloud and n 7 is the number density of 
cloud photons. This rate is to be contrasted to the ordi- 
nary scattering rate of the impinging photon, as derived 

i n J uj 3 uj3. In 



This polarization precession, from an effective 
anisotropic index of refraction, originates in the coherent 
interaction through forward scattering of a single beam 
photon with a large number of cloud photons. In this 
note we develop the theory of another collective interac- 
tion, now between two clouds of photons, also depending 
on coherent forward scattering. This interaction can lead 
to helicity changes when photons of both clouds initially 
all have the same helicity, and to depolarization of one 
cloud when the other cloud is initially unpolarized. The 
rate will now turn out to be of order r po i divided by a 
slowing factor log(iV) where N is the number of photons 
in a region of interaction of linear dimension 1/T p . 

To rederive the Kotkin and Serbo result, and to lay the 
groundwork for the extension, we consider the complete 
set of momentum states {qi} that are occupied in the ini- 
tial state in either cloud (whether singly or multiply oc- 
cupied). We take Lj of and truncate it by keeping the 
parts of the fields that contain only creation and annihila- 
tion operators for this set of momenta. The momentum- 
conserving processes described by this interaction are just 
the forward scattering of beam photons from cloud pho- 
tons since co-moving cloud particles (or beam particles) 
do not scatter from each other in the interaction, l|T]l. 
The result, for the effective "forward" Hamiltonian of the 
system, after substitution of the canonical expressions for 
the electromagnetic fields in terms of creation and anni- 
hilation operators into — Li of Q and performing the 
space integral over a quantization volume, V, is, 



-(11/3) l^l^] 



(2) 



from the cross-section 0, T s = .014 x a 4 m" 
all situations in which uj lj c 
magnitude greater than T s 



where G — 2a 2 /15 m , and where the indices j and m 
extend over the momentum states defined above. 

In (0) the products of photon annihilation and creation 
operators for the beam modes, aj, ajf, and for the cloud 
<< m , T p is many orders of modes, b®, tft (where x and y indicate the polarization 

state and j enumerates the set of momenta) have been 
reexpressed in terms of the operators, 
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T p = ( f)tof -{af )taf . (3) 

with the parallel set of definitions for the cloud opera- 
tors, taking a — > 6, r — > f. The operators r^ 1 ^ 3 ) /2, sup- 
plemented by an operator r^ 2 ^/2, which will not explic- 
itly enter below, obey angular momentum commutation 
rules, as do the operators £/2. Since iJf or only connects 
states of identical unperturbed energies, we did not in- 
clude a contribution from an H in 0. 2 To follow the 
polarization of a single beam photon of energy ui pass- 
ing through the cloud we can then write the Hcisenberg 
equations for f(t) coming from J2J) as, 

4r(t) = -2GUJV- 1 fit) x (Z(t) - vZ^) . (4) 
at 

where we have defined Z — J2 m w mCm > and introduced 
a vector v, defined as a unit vector in the 2 direction in 
the internal space 3 . For the case of an isolated beam 
photon interacting collectively with the cloud photons it 
is fairly clear that we can replace the cloud operators, 
^(!).(3) by their expectation values, since the back re- 
action from beam interactions affects the cloud almost 
not at all. If the cloud polarization is at an angle to 
the x axis and the cloud energies are reasonably nar- 
rowly clustered around an energy uj c we have (Z^) /V = 
LJcHj cos(26>), (Z^^/V — Lu c n 7 sin(20). Since the eqs.Q 
are now linear in the operators for the beam particle, 
they hold for expectation values. Taking the initial con- 
dition (t( 3 )(0)} = 1, (tW<( 2 )(0)} = 0, for an initial beam 
polarization in the x direction, and solving we ob- 
tain the the x, x component of the polarization density 
matrix, P x = 1/2 + (t^ (t)) /2, 

P x = 1 - \ sin 2 (20) [1 - cos(IV )] , (5) 

where T p = 2Gujuj c n-y. 

EqEl recaptures the effects noted by Kotkin and Serbo 
, and we refer the reader to their articles for more dis- 
cussion as to the possibilities of observations. To make 
one comparison to laboratory parameters, we define os- 
cillation length as A = (Tp) -1 and express in ordinary 
units, 

x 1 k lf W£crit\ 2 /lMeV\ 

a = 1.5x10 9 (— ) (— ) cm < 



2 We have made the transition from interaction Lagrangian to the 
interaction Hamiltonian in without regard to the fact that 
the original Lj contains time derivatives of the original canoni- 
cal coordinates, A. This is consistent if we take u>j and u> m as 
simple parameters in and do not translate them back into 
time derivatives when we go to a Heisenberg picture. 

3 One should not confuse the three dimensional space that we have 
introduced using a spinorial representation for the polarization 
vectors with the three dimensional configuration space; we label 
vector components in the former with (1,2,3) and in the latter 
with (x,y,z) 



where B cr it = m 2 c 3 /eh and E is the rms electric field of 
the cloud. In the u>\ = 2.35eV laser used in the experi- 
ment reported in ref.Q the field strength was E/E CT n fh 
1.5 x 10 -6 . In this case taking hu — lOOMeV leads to 
an oscillation length of « 3 cm (The pulse length for 
this laser is a fraction of a millimeter; the free path for 
ordinary photon scattering from the cloud under these 
conditions is of the order of 10 9 cm.) 

We further note that this photon-cloud interaction pro- 
duces no effect, on the short time scale, if the inital po- 
larizations are perpendicular, and we note that if the 
cloud is unpolarized then there is no depolarization of 
the beam. Turning to the case of two colliding clouds, for 
which neither of these conclusions will hold, we assume 
for simplicity that photon densities in the two colliding 
groups are equal. Now we need to take the variables Q on 
the R.H.S. of J5J) as well as the variables f to be dynamic 
variables, rather than taking their expectation values in 
the initial state. 

This calculation is simplest in a helicity basis, however. 
The forward interaction, Hf m , gives a matrix element for 
the transition in which a state of a positive helicity pho- 
ton from one bath and a positive helicity photon from the 
other bath makes a transition to a state with two pho- 
tons of negative helicities. We can easily express Hf OI of 
P]l now in terms of operators £, f) which act in the two 
dimensional helicity spaces of the respective clouds, des- 
ignated respectively as the "up" cloud and the "down" 
cloud. The components l;^ and -qf^ measure the spins in 
the ±z direction for the photons in the respective clouds, 
thus the negative of the helicity in the case of the down- 
moving photon. We choose both clouds to be essentially 
monoenergetic, with energies uj and lo c for the respec- 
tive up-moving and down-moving clouds; then we can 
express the "forward" Hamiltonian in terms of the collec- 
tive coordinates, = J^i^i^i = ^2iVi\ where 
£(+) = (£(i)+i£(2))/2etc. 

By direct transformation of J2J) we obtain, 

H ioI = GwWeV-^+ty - * + 2£(-y+) 

-(11/3) lW lW]. (7) 

Now we pose the question of what happens beginning 
with an initial state in which all N up-moving photons 
have spin +1 and all N down- moving photons have spin 
— 1 in the z direction. We can proceed, as in the earlier 
case, by writing the equations of motion, 

±j(+\t) = -2iGu;uj c V- i e ) (t)ri i+) (t), 
at 

j f ^(t) = -2iGa;a; c F- 1 [e (+) (i)'7 ( - ) (t)-^ ( - ) (*)^ (+) (*)] ) 

j t &\t) = 2iGwu, c V-H i3) (t)r 1 ^(t) , (8) 

plus the three equations in which r and £ in © are inter- 
changed. In the calculation leading to (J5J we proceeded 
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to a soluble problem by taking a factorized ansatz that is 
equivalent, in our present problem, to the replacement, 

« (3) (*)'7 (+) (*)> = « (3) (*)><'? (+) (*)> • (9) 

But for the initial state that we are now considering all of 
the mixing operators with ± superscript have expectation 
value zero, and it is clear that there would be no evolution 
in time at all, were the factorization ansatz valid. We 
proceed instead to a calculation equivalent to solving the 
full coupled operator equations. 

The total z component angular momentum in the new 
internal space in which hclicity is the basis, measured by 
(£( 3 ) + ?/ 3 ))/2 , is conserved. Thinking of the system as 
an assemblage of spins with an upper tier of N spins all 
initially pointed up and a lower tier all initially pointed 
down, we enumerate the states that are connected to the 
initial state (and to each other) by the Hamiltonian of 
|7|. Any number of the N spins in the upper tier, all 
initially up, may be flipped, leading to N + l possibilities 
for the magnetic quantum number of the this tier. The 
operators £-£/4 and ff-ff/4 are separately conserved, each 
with eigenvalue (N/2 + l)N/2. Therefore for each value 
of (^ 3 ^ /2) in our set, there is a single upper tier configu- 
ration that enters, and a single lower tier configuration as 
well. We index the states by the number of flips plus one, 
i, where i takes on the values 1, 2..JV + 1. We express 
the operator products that occur in the Hamiltonian in 
this basis, 

(i\£-r)+\i -l) = (N-i+ l)(i); i = 1, ..N + 1 , 
(t + l|f+»7-|i) = (N — i+ 2)(i - 1); i = 1, ..N + 1, (10) 

which come directly from the standard angular momen- 
tum matrices. We solve numerically for a n + 1 compo- 
nent wave function \&(t), using the Hamiltonian J7J with 
the substitution l|10|) and the initial condition 4"i(0) = 
5i t x, and then calculate the measure of average helicity of 
the upper tier, 

N N+l 

R(t) = N- 1 = J2 l*,W| 2 (^ ~ 2* + 2JJV- 1 . 

i—i z—1 

(ii) 

We perform these calculations for a series of values of 
N, and show the results as a function of scaled time, 
s = T p t = 2GNu c ujt/V. Fig.l displays results for values 
of N ranging from 8 to 512, equally spaced in log(iV). 

The data shown in the figure clearly suggests a charac- 
teristic time of order T" 1 log(iV) for a complete turnover 
of the spins. We can gain a heuristic understanding of 
these results. Instead of the set of operators £, fj we in- 
troduce the bilinear forms, 

x = i^ + ^rj^ ; u = i^'rf + ' ; y = r)^r](~\ 

z = £(")£(+) ; W = S<.3) (12) 

Writing the Heisenberg equations of motion for these op- 
erators by taking commutators with Hf m in the form Q 
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FIG. 1: The function R(s/T p ) of 1111 . the mean helicity of the 
up-moving cloud, for values of N = 8, 16, 32, 64, 128, 256, 512 
as determined from solutions of ®, plotted against the di- 
mensionless scaled time, s The curves for higher values of N 
lie progressively farther to the right. Equal spacings of the 
curves indicates a transition time increasing as log(iV). 

and making the further substitution ?/ 3 ) = — £( 3 ) we ob- 
tain the closed set, 

Nx = w(z + y) — w 2 ; ii = —x ; Ny = wx — uw, 

N z = xw — wu ; Nw = —x + u . (13) 

where the derivatives are with respect to the scaled time 
s. Treating these equations as c-number equations with 
the initial conditions x = y = z = u = and w = N 
allows us to write a single equation for w = w/N, 

d 2 „ , o x 2w 2 . 
— w = 2w(w 2 -l) + — . (14) 

The initial condition is now w(0) = 1, w'(0) = 0. In fig. 2 
we plot solutions of (jl4[) and compare with the numerical 
solutions to the complete equations. The fit is good for 
values of R > .6 . We also see that the in the case of 
solutions to (|14fl the equal spacing continues to values of 
N w 10 4 , leaving little doubt of the logarithmic depen- 
dence. It is possible to understand this limit analytically 
from 1)14(1. capitalizing on the fact that when N — ► oo, the 
solution is the familiar kink solution in a A<fi theory in 
one dimension, w = tanh[(t — to)/2)], then showing that 
for large N, in the time region in question, the w 2 /N 
term in l|14JI can be dropped in favor of changing the 
initial value of w to 1 — 2/N, this in turn determining 
< = log(2iV). 



4 This is exactly the arbitrary procedure that we disparaged above 
for the case of the equations for the operators £ and fj. One 
difference, and perhaps the key to the agreement with results of 
the complete solutions, is that the operators 1121 keep us within 
the TV + 1 dimensional subspace of states defined above, while 
the operators £ and fj do not. 
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FIG. 2: The function R(s/T p ) as determined from the 
solution for the heuristic equation 1141 for values N = 
8, 32, 128, 512, 2048, 8192 (solid lines). The dashed lines show 
the solutions of the complete equations of motion 1111 . as 
plotted in fig.l for the first four values of N. 



To summarize briefly: in many-body systems in which 
every partice of set A interacts with every particle of 
set B, evolution times for macroscopic properties may 
be much faster than one would have predicted based on 
cross-sections, even in the absence of initial phase rela- 
tions among the components that one might have antic- 
ipated were necessary for such behavior. In the photon- 
photon system the effect is an extension of the known 
index-of-refraction effects of photon polarization treated 
in refs. In the detailed example treated, there is total 
oscillation back and forth between all positive helicitics 
and all negative helicities in both clouds. 

The case in which one cloud with 100% polarization in 
helicity collides with an unpolarized cloud is somewhat 
more complex. Here we predict partial depolarization of 
the polarized cloud. From Q we see that photons in 
the target cloud with the opposite helicity to those of 
the beam cloud are effectively sterile. Therefore we can 
discuss the polarization changes of the beam cloud in 
a manner similar to that of the calculation given above. 
There remains an order T p / log(iV) rate of depolarization 
after averaging over the configurations of polarization of 
the individual photons in the target cloud. This is in 
contrast to case of a single photon interacting with the 
cloud discussed at the beginning of this paper, where 
depolarization takes place on the time scale 1 /T s . 

Our calculation was for an idealized system of plane 
wave modes in a box, with (implicit) periodic boundary 
conditions. Does it apply to realizable systems in which 
the two clouds are in contact for a time of order (box 
size/c)? It is clearly required that the characteristic time 
for transformation be shorter than this contact time, a 
criterion that is easily checked in any given situation. It 
is harder to answer the question, "Can the laboratory 
photons in the two beams really sustain a coherent inter- 



action over the whole of the macroscopic region (of order 
of a cm., in the numerical example mentioned above, but 
now multiplied by a logarithm of the order of 100) for 
our process to unfold?" We do not know how to address 
the exact quantification of this question, although we be- 
lieve that the answer is "yes" for the case of the beams 
from lasers and from synchrotrons. Another question is 
that of the role of all of the modes that we have left 
out in using the truncation that produced the "forward" 
Hamiltonian (J2J). We anticipate that over the time-scale 
l/Tp these modes create junk that does not add up to 
anything macroscopically, due to phase oscillations, as 
indeed they must in our preliminary beam-cloud calcula- 
tion. In any case we believe that the "speed-up" through 
the many body interactions that we have described here 
is interesting enough to warrant serious attention to some 
of the harder questions that arise. 

Finally we note the close similarity between the issues 
discussed in this note and in refs. @> @> which discussed 
the possibility of speeded-up flavor transformations of 
colliding neutrino clouds. Although the equations are 
quite similar, a critical difference is a term proportional 
to £( 3 V 3 > on the RHS of the analogue of J7J in the neu- 
trino case. This term destroys the speed-up process in 
the simple model, with just two tiers of states, in which 
all the couplings between the upper tier and lower tier 
neutrino states (in our N-spin terminology) are equal to 
each other. In more realistic (and complex) situations it 
is possible that there would be speeded evolution, how- 
ever. 
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